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Abstract 

We discuss some phenomenological issues of the effective quark- 
lepton operators emerging from the bihnear lepton-Higgs couphngs in 
the superpotential and in the soft supersymmetry (SUSY) breaking sec- 
tor of the supersymmetric models without R-parity. The contribution of 
these operators to the neutrinoless double beta decay (Oz^/3/3) is derived. 

The corresponding nuclear matrix elements are calculated within the 
renormalized quasiparticle random phase approximation, which includes 
the Pauli effect of fermion pairs and does not collapse for the physical 
values of the nuclear force strength. 

On this basis we extract from the experimental data new stringent 
limits on the 1st generation mass parameter characterized the lepton- 
Higgs bilinear coupling and on the electron sneutrino vacuum expecta- 
tion value. 



1 Introduction 

In the standard model (SM) of the electro-weak interactions the baryon B and 
lepton L numbers conservation is protected to all orders of perturbation theory 
by an accidental Uib x Uil symmetry existing at the level of renormalizable 
operators. In the minimal supersymmetric (SUSY) extension of the standard 
model (MSSM) this symmetry is absent and the L and B violating processes 
are not forbidden. A conventional way of eliminating the phenomenologically 
dangerous L,B-violation in this case exploits a discrete symmetry known as 

^ On leave of absence from the Joint Institute for Nuclear Research, Dubna, Russia. 
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i?-parity [0], which is imposed on the modeL This is a multiphcative 
Z2 symmetry defined as i?p = (^— i)3^+'f^+25^ where S", B and L are the spin, 
the baryon and the lepton quantum numbers. R-parity conservation has a 
distinctive phenomenology. It prevents lepton and baryon number violating 
processes, the superpartners are produced in associated production and the 
lightest SUSY particle is stable. The latter leads to the celebrated missing £^7^ 
signature of the SUSY event in high energy detector and renders a cold dark 
matter particle candidate. Although desirable for many reasons the R-parity 
conservation has no well motivated theoretical grounds. 

On the other hand relaxing the R-parity conservation we may get a new 
insight into the long standing problems of particle physics, in particular, to 
the neutrino mass problem. Remarkable, that in this framework neutrino can 
acquire the tree level supersymmetric mass via the mixing with the gauginos 
and higgsinos at the weak-scale [^, [0-[[7[. This mechanism does not involve 
the physics at the large energy scales Mint ~ C(lO^^GeV) in contrast to the 
see-saw mechanism but relates the neutrino mass to the weak-scale physics 
accessible for the experimental searches. 

The R-parity can be broken ) either explicitly or spontaneously 
The first option allows one to establish the most general phenomenologi- 
cal consequences of R-parity violation while a predictive power in this case is 
rather weak due to the large number of free parameters. Spontaneous realiza- 
tion of Tjiy SUSY is much more predictive scheme leading to many interesting 
phenomenological consequences However, it represents a particular model 
of the R-parity violation. At present it is an open question which underlying 
high-energy scale physics stands behind the R-parity, protecting or violating 
it at the weak scale. 

Many aspects of the ^pSUSY models in high and low energy processes had 
been investigated in the literature [|16]-[^ . 



Recently, a growing interest to the supersymmetric models without R- 
parity was stimulated by the exciting news from the HERA experiments, re- 
ported the anomaly in deep inelastic e"'"p-scattering |jl5[ which can be elegantly 
explained within these theoretical framework in terms of the lepton number 
violating interactions. 

Since the lepton number is not conserved without R-parity some low-energy 
exotic processes become possible within the MSSM. Among them the neu- 
trinoless nuclear double beta decay (Oz//5/5) is known to be very sensitive to the 



certain ^pinter act ions |]T8|. Provided an unprecedented accuracy of the modern 



0z//9/3-decay experiments |21| this allows one to establish stringent constraints 
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on the ^pSUSY ]T6[-|20 



In the present paper we consider the imphcations of the bihnear lepton- 
Higgs Ijlp terms on Oz//3/3-decay. In the general case of the explicitly broken 
R-parity these terms are present in the superpotential and in the soft SUSY 
breaking potential. Previously the main attention was paid to the phenomenol- 
ogy of the trilinear Yukawa couplings. It was widely believed that the bilin- 
ear terms can be rotated away by a proper field redefinition. However, it is 
not the case in the presence of the soft SUSY breaking interactions 0, 0. It 
was realized that the bilinear violation, generically leading to the non-zero 
vacuum expectation values (VEV) of the sneutrino fields and to the lepton- 
gaugino-higgsino and slepton-Higgs mixing, provides a number of interesting 



phenomenological issues in 



In particular, this mixing generates the new effective lepton number vio- 
lation operators which contribute to the nuclear Oz//?/3-decay. In what follows 
we derive these operators and analyze their net effect in the presence of the 
nuclear media. 

The paper is organized as follows. Basic ingredients of the Iflp MSSM with 
the general setting of the explicit R-parity violation are shortly described in 
Section 2. In Section 3 we discuss the bilinear ^p mechanism of the nuclear 
Oi//3/3-decay. Here we analyze all the tree-level -^pMSSM contributions to the 
Oi//3/?-decay amplitude. We start with the quark level and derive the corre- 
sponding low energy effective Lagrangian. In Section 4 we take into account 
the effect of nuclear structure and derive the corresponding nuclear matrix el- 
ements. Then we calculate their values within the renormalized Quasiparticle 
Random Phase Approximation (pn-RQRPA) [^]. The pn-RQRPA is an ex- 
tension of the pn-QRPA by taking into account the effects of the Pauli principle 
for the fermion pairs. In this approach the sensitivity of the nuclear matrix 
elements to the details of the nuclear Hamiltonian is reduced considerably. 
Using experimental lower bound on the ^^Ge half-life we extract in Section 5 
stringent constraints on the 1st generation lepton-Higgs mixing mass parame- 
ter and on the electron sneutrino VEV. We close our discussion with the short 
comments on some implications of these constraints for the other experiments. 
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2 Minimal SUSY model with R-parity viola- 
tion. 



In order to set up our notations let us briefly recapitulate the main ingredients 
of the minimal SUSY standard model (MSSM) with explicit R-parity violation 
(^pMSSM). 

The Rp violation is introduced into the theory through the superpotential 
and soft SUSY breaking sector. 

For the minimal MSSM field contents the most general gauge invariant 
form of the renormalizable superpotential reads 

W = Wr^ + W^^. (1) 

The Rp conserving part has the standard MSSM form 

y^R, = HlHiLE' + HdHiQD' + huH2QU' + fiHiH2. (2) 

Here L, Q stand for lepton and quark doublet left-handed superfields while 
E'^, f/^, D'^ for lepton and up, down quark singlet superfields; Hi and H2 
are the Higgs doublet superfields with a weak hypercharge Y = —1, +1, 
respectively. Summation over the generations is implied. 

The Rp violating part of the superpotential (0) can be written as 0, § 

= XijkULjEl + X'iji^LiQjDl + fiiLjH2 + X'-j^UiDjDl, (3) 

The coupling constants A (A") are antisymmetric in the first (last) two indices. 
The first two terms violate lepton number while the last one violates baryon 
number conservation. 

Another source of the R-parity violation is the soft supersymmetry breaking 
part of the scalar potential. It contains the -terms 



v;;^' = k,kL.L,Ei + ~k,,UqM + K.^-upi + ~^^UH2+ (4) 

+m\HLiH\ + H.c. 

The simultaneous presence of lepton and baryon number violating terms in 
Eqs. (0), (|) (unless the couplings are very small) would cause unsuppressed 
proton decay. Therefore, either the lepton or the baryon number violating 
couplings can be present. There may exist in the theory an underlying discrete 
symmetry such as the B-parity 0, ||23| which forbids dangerous combinations 
of these couplings. Henceforth we simply set A" = A" = 0. 
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The remaining R-parity conserving part of the soft SUSY breaking sector 
includes the scalar field interactions 




(5) 



i=scalars 



-huAuH2QU' - fiBHiH2 + H.c. 



and the "soft" gaugino mass terms 



- \M-^BB + M^W^W^ + M^g^-g"] - H.c. 



(6) 



As usual, M^^2,i denote the masses of the S'[/(3) x SU{2) x f/(l) gauginos 
g, W, B while mj stand for the masses of the scalar fields. The gluino g soft 
mass M3 coincides in this framework with its physical mass denoted hereafter 
as Trig = M3. Al, Ad, Ajj and B in Eq. are trilinear and bilinear "soft" 
supersymmetry breaking parameters. All these quantities are free SUSY model 
parameters which due to the renormalization effect depend on the energy scale. 

In this paper we assume for simplicity the universal gaugino soft masses at 
the grand unification scale Mqut- At the weak scale this leads to the following 
relations 



An impact of the R-parity violation on the low energy phenomenology 
is twofold. First, it leads the lepton number (LNV) and lepton flavor (LFV) 
violating interactions directly from the trilinear terms in IV^p. Second, bilinear 
terms in Wfi^ and in V^"'^* generate the non-zero vacuum expectation value for 
the sneutrino fields (z/j) 7^ and cause neutrino-neutralino as well as electron- 
chargino mixing. The mixing brings in the new LNV and LFV interactions 
in the physical mass eigenstate basis. Below we will specify those interactions 
which are relevant for the Oz//3/3-decay. 

The trilinear terms of the R-parity breaking part of the superpotential W^i^ 
lead to the following AL = 1 lepton-quark operators 

-^^A = >^ijk[i>iLekPLej + ejiekPLi^i + ekRejPRi^i - {i ^ j)] + (8) 



Ml = (5/3) tan 9^M2, M2 ^ 0.3M 



(7) 



+ KjA^iLdkBidi + djidkBh^i + dkRdjPRU- - endtPLUj 
- UjLdkPiGi - dkRUjPRef\ + H.c. 



Here, as usual Pl,r = (1 =F 75)/2. 
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The presence of the bihnear terms in the Eqs. (0),(§) leads to the terms in 
the scalar potential linear in the sneutrino fields z/j. As a result, at the mini- 
mum of the potential (i>j) 7^ 0. Thus, the MSSM vertices Zvv and Wev create 
the gaugino-lepton mixing mass terms Zv{y) e{v) (with W being wino 
and zino fields). Combining this terms with the lepton-higgsino HiLiHi mixing 
from the superpotential Eq. we end up with 7x7 neutral fermion and 5x5 
charged fermion mass matrices (see Appendix A). The mass eigenstate fields 
can be written in the form 

^(0). = S,, ^'(0),, ^(±)^ = (9) 
with the weak eigenstate fields in two component notation 

vl/'j) = {ui, -iX', -^3, Hi Hi), (10) 
^'f) = (eZ, /iZ, rl, iff), (11) 

^(+) = (el, rt. -iK. Ht). (12) 

Here z/j are the neutrino fields. A' and A3, A_ are the Uiy and SU^l gauginos, 
respectively while higgsinos are denoted as Hi2iHf2- The mixing matrices S 
and diagonalize the neutralino-neutrino and the char gino- charged lepton 
mass matrices respectively. The lightest mass eigenstates are identified with 
the physical neutrinos and the charged leptons. Remarkable, that as a result 
of the minimal field content and the gauge invariance the neutral fermion mass 
matrix M.q ( |A.4|) before diagonalization has such a texture that its first three 



rows and the last one are linearly dependent and, as a result, two neutrino mass 
eigenstates are degenerate massless states. The third neutrino state acquires 
the tree level mass which approximate form is (see Appendix B) 

It is natural to identify the massive neutrino state with the tau neutrino Vr 
while the two massless states with the Pe and v^. The Ve — niass degeneracy 
is lifted by the 1-loop corrections as well as by the non-renormalizable terms 
in the superpotential giving to z/g,^ the small non-equal masses [^. As to the 
tau neutrino mass in Eq. ([13| ) it is subject to the experimental constraint 
f^ur ^ 23MeV |2^. Assuming no cancellation in Eq. (|13]) this leads to the 
upper bounds 

/ii < 15GeV, (z>,) < 7GeV. (14) 
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at the typical sample values of the MSSM parameters /i ~ M2 ~ M^. Of 
course, these bounds are only indicative and may essentially vary from point 
to point in the MSSM parameter space. 

The rrij,^ constraints can be evaded assuming an approximate alignment 
between two vectors = and 6, = ((t'i), {Hf)) which leads to the can- 

cellation in Eq. (|13D since . This might be guaranteed 

by a special global symmetry or by some dynamical reasons @]. 

Rotating the MSSM Lagrangian to the mass eigenstate basis one obtains 
the RPM generated lepton number violating interactions which bring many 
interesting implications for the low and high energy phenomenology. Below 
we are studying they contribution to the Oz//5/5-decay. 



3 LH-induced Oupp-decay. Quark level transi- 
tions. 

We have analyzed all the possible tree level contributions to the Oi//3/3-decay 
amplitude which include the RPM interactions and the superpotential A, A' 
couplings from Eq. (P). The leading diagrams are presented in the Fig. 1. 
The diagrams in Fig.l(a,b) incorporate only the RPM generated vertices, and 
in Fig. l(c,d) these vertices are accompanied by one A' type vertex (on the top 
of the diagrams). The diagram in Fig. 1(a) has in the intermediate state either 
neutrinos or neutralinos and two W-bosons while the diagrams in Fig. l(c,d) 
neutrinos, squarks/selectron and one W-boson. The diagram in Fig. 1(b) is 
mediated by the gluino and double squark exchange. The diagram Fig. 1(a) 
with the neutrino exchange is the conventional Majorana neutrino contribu- 
tion to the Oz//3/3-decay. Recall that in the ^pMSSM with the bilinear R-parity 
violation the neutrino masses and mixing angles are derived at the tree level 
in terms of //j, (z/j) and the MSSM parameters (see Appendix A). Therefore, 
this contribution inherently pertains to this model. We did not include in this 
list those diagrams which do not contain RPM vertices. These diagrams con- 



structed of the A, A' couplings were previously analyzed in Refs. |]T6|]-|2y]. All 
the other diagrams in this order of perturbation theory have extra suppression 
factors and, therefore, can be neglected. The suppression factors originate from 
the smallness of neutrino mass, when it appears in a positive power, from the 
1st generation left-right sfermion mixing proportional to mu,d,e/MsusY and/or 
from the fermion-sfermion-higgsino couplings proportional to mu,d,e/Mw with 
^« d e being the u,d quark and the electron masses respectively while Msusy 
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denotes the typical SUSY breaking mass scale. 

Now let us specify those RPM generated operators which are encountered 
in the diagrams in Fig. 1. They are 

i:-Ln = -^KnW-erPLXn + (15) 

+V2g2 {Pp.PRd dl + PlvkPRu'' + PliVkPRe" e,^ + CuPrc' d,) + H.c. 

The subscripts k, i denote generations. 

The first term is generated from the standard model W — e — v and the 
MSSM W — ~X interaction terms while the rest originates from the MSSM 
neutralino (chargino)-fermion-sfermion interactions X ~ Q ~ X^ ~ Q ~ Q (fo^' 
the MSSM Lagrangian see 

Note that the trilinear fermion-sfermion couplings in Clh are not present 
among the superpotential trilinear A, A' terms in Eq. (|^). 

The coefficients in Eqs. (|T5[) depend on the mixing matrix elements intro- 
duced in Eq. (H): 

3 

/^n = X! ^ll'='n+3,3 + A;^4.=:,„+3^5 + (16) 
k=l 

= — 77 (fan6',^^fc4 + S^fcs) ! 
o 

d 1 1 

i^k = -g tan6'yi/SA:4, C = -—^A^4^. 

In what follows, for the derivation of the constraints on parameters (z/j),/ij 
characterizing the bilinear ^pwe employ the approximate analytical diagonal- 
ization method of the Ref. [|1^]. It allows one to represent the mixing matrices 
in a convenient analytic form and express the dependence of the coefficients 
in Eqs. (|1^) on the afore-mentioned parameters explicitly. In the leading 
order in small parameters fii/Mz, {vi) /Mz we obtain 

= iluKk-^illK2-il2K3, (17) 

L 



f^k = - V/;)* (tan^^eii + 3^*2 



Pt = -tan^^4^)*e;i, C=^en. 
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The notations used in these formulas are explained in Appendix B. 

The MSSM gluino-quark-squark vertex in the diagram Fig. 1(b) is de- 
scribed by the Lagrangian term 

C~, = -V2g,^ {rPR~9^''^l - rPL-g^^Vn) + h.c, (18) 

Here A*-"-* are 3x3 Gell-Mann matrices (a = 1, 8). Superscripts a, (3 denote 
the color indices. 

The diagrams in Fig.l describe the Ijip SUSY induced quark transitions 
which proceed in the nuclear media and trigger the nuclear Oi//3/5-decay. Our 
goal is to derive the corresponding half-life for a certain isotope assuming for 
simplicity that there is no other contributions to this nuclear process. In order 
to apply the standard approach |]18[, based on the non-relativistic im- 



pulse approximation one has to derive first the effective low energy Lagrangian 
describing the basic Oz//3/3-quark transition dd — > uuee in terms of the color 
singlet quark charged currents which can be embedded then into the corre- 
sponding hadronic (nucleon or pi-meson) currents inside a nucleus. One has 
also to separate the short and long ranged parts of the quark level transition 
operators since they are treated within this approach in different ways. It 
is understood that the short ranged parts involve only heavy particles in the 
intermediate states {XtX^ ^^A^^) while the long distance ones include the 
neutrino exchange. 

Integrating out the heavy fields from the diagrams in Fig. 1 and carrying 
out Fierz reshuffling we obtain the desired effective Lagrangian which allows 
one to reproduce the low energy contribution of these diagrams in the first or 
in the second order of perturbation theory. It takes the form 

- V2GfK,, ■ # ) ■ {ukPne') J + GFV2{n''PLi^k)vlk^ J,. 

Here we introduced the color singlet quark currents 

J = u''PRda, J^'"" = v!'a^''PRd^, J^' = MV^Lt^a, (20) 

The effective parameters r] accumulating the dependence on the initial 
SUSY parameters are defined as 



(eP^e^) - (19) 
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(ki) 



Vx = 7^ 2 dii dii 



Here rrig, rriq and m^. are the gluino, squark and neutralino masses. 

In the MSSM we have for the neutrino mixing matrix element (see 
Appendix B) the following expression 

= Sik cos e - 63k sin e (22) 

with 

sine = -Ai/|A|. (23) 



In Eq. ([T9|) the first and the second terms reproduce the contribution of 
the gluino Fig. 1 (b) and the neutralino Fig. 1(a) exchange graphs in the 
1st order of perturbation theory while the third and the last terms reproduce 
the contribution of the neutrino exchange graphs in Fig. 1(a) and Fig. l(c,d) 
in the 2nd order of perturbation theory. Note, that the last term in Eq.(|l9D 
is the ordinary lepton-number L conserving standard model interaction term. 
Since the 0t//9/9-decay requires the L-violation AL = 2 the neutrinos in Fig. 
1 (a) propagate in the Majorana lepton-number violating mode. In this case 
the source of the L-violation is given by the neutrino Majorana mass term. 
That is why the contribution corresponding to this diagram is proportional 
to the neutrino mass or, more precisely, to the average neutrino mass (m,y) 
defined below. On the contrary, the neutrino exchange diagrams in Fig. 1 
(c,d) are not proportional to {m,y) and survive in the limit rriy = since the 
lepton-number violation AL = 2 is produced by the interaction term vPrc'^ in 
Eq. ( [TP| ) itself. Therefore, the neutrinos in diagrams Fig. 1 (c,d) propagate in 
the L-conserving Dirac mode. 

So far we concentrated on the Oz//3/3-transitions at the quark level described 
by the effective Lagrangian ([191) . The aim of this paper is the calculation of the 
amplitude for the nuclear Oi^/5/3-decay taking into account nuclear structure. 

The next section deals with the derivation of the amplitude for the nuclear 
Oi//5/5-decay triggered by the quark transitions in Fig.l. 



4 Nuclear Oz//3/3-decay. 

Let us write down the following formal expression for the amplitude Rquib/s of 
Oi//3/3-decay 

7^o./3/3 = <iA,Z + 2),2e-\S-l\iA,Z)>= (24) 
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< {A,Z + 2),2e~\Texp[i J d'^xCeffix)]\{A, Z) > 



where the effective Lagrangian Ceff is given by Eq. (pj]). The nuclear structure 
is involved via the initial (A,Z) and the final (A, Z+2) nuclear states having 
the same atomic weight A, but different electric charges Z and Z+2. The 
standard framework for the calculation of this nuclear matrix element is the 
non-relativistic impulse approximation (NRIA) 

It is straightforward to derive the following formula for the amplitude of 
the 0+ 0+ transition with two outgoing S-wave electrons 



7^, 



01//3/3 



(0+^0+) =Co./le(l + 75)e^x 



(25) 



rrir, 



< > 



-Mn + 



rrip 



The normalization factor is 

Cou = {Gl2m,)/{8V27rR). 



(26) 



Here, rrie and Ja ~ 1.261 are the electron mass and the nucleon axial coupling, 
R = tqA^^^ is the nuclear radius {tq = 1.1/m). 

The last term is the conventional Majorana neutrino mass contribution 
proportional to the average neutrino mass. In the MSSM we have 



2 _ 2 glM2 ^2 



3 DetM, 



(27) 
and 



Here we neglected the small loop induced neutrino masses rriu^ ~ m,^^ 
used Eq. (H). 

Let us specify the nuclear matrix elements involved in the formula for the 
Oi//5/3-decay amplitude (|^). They are 



Mn 



m, 



m 



m, 



ap 
rripR 



1 



-M 



GT,X 



+ M 



(28) 



\mpj me^yfA/ 
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Here and me stand for the proton and electron masses, fy ~ 1-0 is the 
vector nucleon constants, = 0.85GeV to be defined below. The coefficient 
ap = 1.75 is related to the nucleon matrix element of the pseudoscalar current. 



Its numerical value calculated in the quark bag model we take from Ref. [18 



The partial nuclear matrix elements in the closure approximation we write 
down in the form 



MF,i 
MGT,i 
MT,i 



(0;iE^aV^.M ( — t Oab |0+), 



(29) 



(0;iE^aV^(^a^>) 



-)' 

Tab J 



Sab \^t)-> 



where i = g, X, N, v. The exponent takes the values 5i = {1, 0, 1, 0}. We use 
the shorthand notations 

J^i = {0,0, FN{xA),h+{rab)}, % = {F2{x^),hT'{rab),0,0}, (30) 
Gi = {Fi{x^),hR{rab),FN{xA),h+{rab)}, 

for the following form factor functions and neutrino potentials 



Fi{x) 
F2{x) 

h+{rab) 
hnirab) 
hr'irab) 



i^3 + 3x + a;2 



a 



e-^ Fn{x) = ^(3 + 3x + x^)e 



+ a2"(x + l) 



x^ 



2 

TT 

2 i?2 



/•OO 

R\ dq- g$2(q2 
Jo 

TT m,, JO 

Jo 



ab 



2 i?2 



q + A ' 

Jojqrgb) 
q + A ' 



(31) 



TT 3m, 



q + A 



with g = |q| being an absolute value of the 3-momentum transferred between 
the decaying nucleons. a^'" = —4.4 ■ 10~^ and a^'^ = 0.2 are the pion structure 
coefficients introduced and calculated in Ref. PU|. A ^ lOMeV is the aver- 
age excitation energy of the intermediate nuclear state. The spherical Bessel 
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functions are defined in the standard way 



MX) 



smx 



sm X cos X 



J2{x) = -ji{x) - jo{x). (32) 



The nucleon form factor in Eqs. (|3TD takes into account the finite nucleon 
size. In our numerical analysis we employ the conventional dipole parametriza- 
tion 



$(q2) 



1 



(33) 



with ruA = 0.85GeV. We also defined in Eqs. (|3g)-(|2D: 



Sab = 3((r„ ■ f„J ((7^ ■ f„J -a^ ai,, 

fab = (fa fb)j ^ab = I'^abj; ^ab 
Xa = f^A^ab-, = Trij^Tab-, 



aab = a^-a^, 

■ Tab /fab, 



(34) 



where is the coordinate of the "ath" nucleon. 

The following comments on the nuclear matrix element A^g in Eq. 
associated with the gluino graph Fig. 1(b) are in order. As discussed in Ref. 

it consists of the two parts -Mg = + corresponding to the 

gluino graph contribution via the two-nucleon and the pion-exchange modes 
respectively. These two modes arise from the two possibilities of hadronization 
of the 1st term of the effective Lagrangian Cgff in Eq. (|19|). One can place 
the four quark fields present in this term in the two initial neutrons and two 
final protons separately (2N-mode). Then nn —>■ + 2e~-transition is directly 
induced by the underlying quark subprocess dd ^ uu + 2e~. In this case the 
nucleon transition is mediated by the exchange of a heavy particle which is the 
gluino g with the mass rrig > lOOGeV. Therefore, the two decaying neutrons 
are required to come up very closely to each other what is suppressed by the 
nucleon repulsion. Another possibility is to incorporate quarks involved in the 
underlying SUSY transition dd ^ uu + 2e~ not into nucleons but into two 
virtual pions or into one pion as well as into one initial neutron and one final 
proton . Now nn pp + 2e~ transition is mediated by the charged pion- 
exchange between the decaying nucleons (vrN-mode). Since the interaction 
region extends to the distances ~ l/ni^^ this mode is not suppressed by the 
nucleon repulsion. An additional enhancement of the vrN-mode comes from 
the hadronization of the SUSY effective vertex operator U'-y^d ■ U'y^d ■ cPrc^ 
replaced by its hadronic image vr^ ■ cPrc^. The enhancement occurs due to the 
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coincidence of the pseudoscalar quark bilinears wy^d with vr-meson field. As 



is shown in Ref. |20| the vrN-mode absolutely dominates over the 2N-mode. 



Therefore we neglected the subdominant 2N-mode part M.^ in Eq. 



We calculate the nuclear matrix elements within the renormalized Quasi- 
particle Random Phase Approximation (pn-RQRPA) [^. This nuclear struc- 
ture method has been developed from the proton-neutron QRPA approach, 
which has been frequently used in the Oz//5/5-decay calculations. The pn- 
RQRPA is an extension of the pn-QRPA by incorporating the Pauli exclusion 
principle for the fermion pairs. 

The limitation of the conventional pn-QRPA is traced to the quasiboson 
approximation (QBA), which violates the Pauli exclusion principle. In the 
QBA one neglects the terms coming from the commutator of the two bifermion 
operators by replacing the exact expression for this commutator with its ex- 
pectation value in the uncorrelated BCS ground state. In this way the QBA 
imphes the two-quasiparticle operator to be a boson operator. The QBA leads 
to too strong ground state correlations with increasing strength of the residual 
interaction in the particle-particle channel what affects the calculated nuclear 
matrix elements severely. 

To overcome this problem the Pauli exclusion principle has to be incor- 
porated into the formalism p2| in order to limit the number of quasiparticle 
pairs in the correlated ground state. The commutator is not anymore boson 
like, but obtains corrections to its bosonic behavior due to the fermionic con- 
stituents. The pn-RQRPA goes beyond the QBA. The Pauli effect of fermion 
pairs is included in the pn-RQRPA via the renormalized QBA (RQBA) p2| , 
i.e. by calculating the commutator of two bifermion operators in the corre- 
lated RPA ground state. Now it is widely recognized that the QBA is a poor 
approximation and that the pn-RQRPA offers the advantages over pn-QRPA. 
Let us stress that there is no collapse of the pn-RQRPA solution for a physi- 
cal value of the nuclear force and that the nuclear matrix elements have been 
found significantly less sensitive to the increasing strength of particle-particle 
interaction in comparison with QRPA results. Thus, the pn-RQRPA provides 
significantly more reliable treatment of the nuclear many-body problem for the 
description of the Oz//3/3 decay. 

For numerical treatment of the Oi//3/3-decay matrix elements listed in Eqs. 
( |30[ ) within the pn-RQRPA we transform them by using the second quanti- 
zation formalism to the form containing the two-body matrix elements in the 
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Table 1: Nuclear matrix elements for the neutrinoless double beta decay 
™Ge(0+) Se{0+) within the pn-RQRPA. 



■Mgt,n 


Mf,n 


MgT,u 


Mf,u 


Mgt,x 


Mt,x 


0.071 


-0.025 


2.6 


-1.2 


1.20 


0.21 


MGT,g 


MT,g 


Mg 


Mx 


M, 


Mn 


-0.34 


-0.089 


-649 


88 


-3.4 


-132 



relative coordinate. One obtains |E6|: 



< Oi2 > 



E 



{2 J + 1) I 

L Jn' 



jp' J 

^ <p.p'-,J\f{ri2)rtrt0^2f{ri2)\n,n'-J> x 
< 0/ II [Cp,c„/]j II J^'nif >< J^'mflJ^mi >< J^'nii \\ [c+c„]j 



0+ > 



(35) 



Oi2 represents the coordinate and spin dependent part of the two body tran- 
sition operator of the Oi//3/3-decay nuclear matrix elements in Eqs. (ppj). The 
short-range correlations between the two interacting nucleons are taken into 
account by a correlation function 



f{r) 



e-°^'Vl 



br 



with a = 1.1 fm^ and b = 0.68 fml (36) 



The one-body transition densities and other details of the nuclear structure 



model are given in |22, 26 



The calculated nuclear matrix elements for the Oi//3/3-decay of A=76 isotope 
within the pn-RQRPA are presented in Table 1. The considered single-particle 
model space has been the 12-level model space (the full 2— 4^ci; major oscillator 
shells) introduced in Ref . . The nuclear matrix elements listed in the Table 
1 have been obtained for the Qpp = 1.0 where Qpp is introduced to renormalize 



the particle-particle interaction strength of the nuclear Hamiltonian. 

According to our numerical analysis, variations of the nuclear matrix ele- 
ments M.g, A4\, A4n and A4i, do not exceed 15% and 30% respectively within 
the physical region of the nuclear structure parameters. 
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Having all the quantities in the Oi//3/?-decay amplitude Eq. (^) specified we 
are ready to extract the limits on the parameters from the non-observation 
of the Oi//3/3-decay. 

5 Oz//?/?-decay constraints on bilinear R-parity 
violation. 

Starting from the Eq. (^3p we derive the half-life formula 

[T^/f /3(o+ ^ 0+)]-^ = GoilA^.I VP- (3T) 

Here Gqi is the phase space factor tabulated for various isotopes in Ref. P7[| . 
We introduced the dimensionless parameter 

(mj^) rrip 

A = \ -f UN +'ngUJg + A^ll 7]^ ujx, (38) 

where a;, = with i = g,X,N. The first, second and third terms in 

this equation correspond to the contributions of the neutrino, neutralino and 
gluino graphs in Fig. l(a,b). Graphs in Fig. l(c,d) contribute to the last term 
in Eq. (^) . It is worthwhile noticing that at typical randomly sampled values 
of the MSSM parameters M2, fi, tan/3 the neutrino exchange contribution from 
Fig. 1(a) dominates over the other contributions. 

The most stringent experimental lower limit on the Oz//5/5-decay half-life 
has been obtained for ^^Ge pTH 

yOv/3/3-exp^Q+ ^0+) > 1.1 X 10^5 years 90% c.l. (39) 

With the nuclear matrix elements calculated in the previous section this 
lower limit can be cast into the following upper bound 

|^|<1.0-10"^ (40) 

This constraint represents a complex exclusion condition placed by the non- 
observation of the Oi//3/3-decay on the MSSM parameter space. The indi- 
vidual bounds on the bilinear parameters /ii, (ui) of our present concern 
depend on concrete SUSY model settings which fix the values of M2, fj, and 
tan/? in the left hand side of Eq. (^0]). 
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Typical constrains for the 1st generation _^pparameters /ii, (ui) , fiiX^n, 
can be obtained at the typical weak scale values of the MSSM parameters 
M2 = 11 = lOOGeV^ and tan/? = 1. We also assume, as is commonly done in 
the similar cases, the absence of a significant cancellation between the terms 
in the left hand side of Eq. (|40| ) defined in Eq. (^). Thus, we come up with 
the following constraints 

< 470KeV, (41) 
|(z>i)| < 840KeV, (42) 
l/UiA'ml < lOOeV, (43) 
|(^>i)A'nil <55eV (44) 

Recall that in our notations {vi) = (i>e). To our knowledge these stringent 
constraints for the 1st generation ^pparameters were not previously considered 
in the literature except a parenthetic note in Ref. 0. One can find in the 
published papers only those constraints which involve the combinations of the 



1st and 2nd generation bilinear parameters or contain only the 3rd 



generation ones [|T3 



To see how stringent are the obtained constraints we can compare them 
with the following one 

A'm < 1.3-10-^ (45) 



which is known as a most stringent constraint on the R-parity violation pO . 
This constraint was previously obtained from the 0z//3/3-decay by taking into 
account only the superpotential trilinear couplings in Eq. (^. Consider for 
comparison the dimensionless quantity 

,.«„^.t^LM, (46) 

with MsusY ~ lOOGeV being the typical SUSY breaking scale. As follows 
from Eqs. (|l5|)-(p!7|) this dimensionless quantity sets the strength of the RPM 
induced trilinear fermion-sfermion-fermion interaction similarly to the coupling 
Am in Eq. (|]). This makes reasonable the comparison of the constraints 
placed on these couplings by the experiment. From Eqs. (|^ ) -(|42| ) we get an 
estimation 
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This constraint looks more stringent (if such a comparison is legitimate) than 
that for Am in Eq. (H). 

After all we conclude that the R-parity violation within the 1st generation 
is restricted by the Oi//3/3-decay to a very low level. Now this statement holds 
for the generic case of the ^pSUSY including both the superpotential trilinear 
couplings and the bilinear terms in the superpotential as well as in the soft 
supersymmetry breaking sector. 

This conclusion has some immediate phenomenological consequences for 
the other experiments, in particular for the accelerator ones. For instance, 
among the two body decay modes of the neutralinos 

and similar processes open in the presence of the bilinear terms one can 
now safely neglect the modes with electron or v^.. 

We can also generalize the arguments used in the SUSY interpretation 
of the HERA anomaly . It is believed that this anomaly can be explained 
by the s-channel squark exchange q\e q* qje, xqi, X^l'i between the initial 
quark-lepton state and the final state particles. The quark-lepton vertex qqe 
allowed in the ^pSUSY models receives the contributions both from the trilin- 
ear A' couplings and from the trilinear operators induced by the bilinear terms 
via the lepton-gaugino-higgsino mixing. It is a common practice to neglect 
the 1st generation squarks in the above mentioned $,p SUSY explanation of 
the HERA anomaly. The argument is derived from the stringent constraint 



on the 1st generation A'j^^^^ coupling pOf shown in Eq. (^51) . However, it does 



not take into account the effect of the bilinear Ifip operators. Now, having at 
hand the new stringent limit on the 1st generation bilinear R-parity violation 
in Eq. (^TD-(PD we can extend the validity of this argument to a general case 



of R-parity violation considered in the present paper. 



6 Conclusion. 

In summary, we derived the contribution of the bilinear R-parity violating 
terms to the neutrinoless double beta decay. Alone with the analysis of the 



trilinear terms previously made in Refs. |1T6[-||20[ this completes the derivation 
of all possible tree- level contributions to the Oz//3/3-decay within the ^pMSSM. 

From the non-observation of 0z//9/9-decay we obtained new stringent up- 
per limits on the 1st generation R-parity violating parameters such as the 
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lepton-Higgs mixing mass parameter and the vacuum expectation value 
of the electron sneutrino (Dg). Then we discussed some implications of these 
constraints on the other experiments and, in particular, on those which are 
running or planned at accelerators. We conclude that the R-parity violating 
effects within the 1st generation, if exist, are very small and in most cases can 
be neglected in phcnomenological analysis of observable effects. 

A special attention was paid to the effects of the nuclear structure in the 
Oi//3/3-decay. In the framework of the pn-QRPA approach we obtained the 
nuclear matrix elements which are stable with respect to the variation of the 
nuclear model parameters within the physical domain. Thus, wc believe that 
our conclusions concerning the particle physics side of the Oi//3/3-decay do not 
suffer from the nuclear structure uncertainties. 



We are grateful to V.A. Bednyakov for helpful discussions. S.K. would like 
to thank the "Deutsche Forschungsgemeinschaft" for financial support. 



Below we present the mass matrices of the neutral and charged fermion sectors 
for the general case of the bilinear R-parity violation within the MSSM field 
contents. 

1.1 Neutral fermion mass matrix. 

In the two component Weyl basis 
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1 Appendix A 



/T 



(0) - 



(A.l) 
(A.2) 



the mass term of the neutral fermions is 



>C^L = -2^(0)^0^(0) + H.c, 



(A.3) 



The 7x7 mass matrix has the distinct see-saw structure 




(A.4) 
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with 3x4 matrix 



m 



I —MzSwC/3Ui MzCwCj3Ui —jJLi 
-MzSwCf3U2 MzCwCpU2 -/i2 
V -MzSwCf3U3 MzCwCpUs 



(A.5) 



originating from the ^pbihnear terms in the superpotential and the soft SUSY 
breaking sector. 

In Eq. ( |A.4| ) is the usual 4x4 the MSSM neutrahno mass matrix in 
the basis {—iX', —iXs, Hi, H2} 



/Ml -MzSwCf3 Mzswsp \ 

M2 Mzcwc,3 -Mzcwsp 

-MzSwCf^ Mzcwcp -/i 

V Mzswsp -MzCwS/B -fi / 



(A.6) 



Here Ui = (ui) / {Hi) and tan/5 = {H2) / {Hi) and sw = cos9yy, cw = cos^^^, = 
sin/3, C/3 = cos/3.. 

In the mass eigenstate basis defined as 

^(0). = S,,^'(o),, (A.7) 

the 7x7 neutral fermion mass matrix TWq in Eq- ( |A.4| ) becomes diagonal 

E*Mor.^ = Diag{m^^,m^J, (A.8) 

where rrii,^ and m^. are the physical neutrino and neutralino masses. For the 
considered case of the tree level mass matrix the only one neutrino has a non- 



zero mass m 



m 



^2 



0, niy^ ^ 0. 



1.2 Charged fermion mass matrix. 

The mass term of the charged fermion sector has the following form 

^Ll. = -^'f )A<±^'(+) + H.c. (A.9) 
in the two component Weyl spinor basis 



(e^, l^L^ ^1 



(A.IO) 
(A.ll) 



20 



The 5x5 charged fermion mass matrix is 

m(') E 

where M^± is the MSSM chargino mass matrix 



-M. = , (A.12) 



The sub-matrices E and E' lead to the chargino-lepton mixing. They are 
defined as 

/ \/2MzCwCpUi /i 

E= y/2MzCwCpU2 /i2 |, (A.14) 
V V2M^Ch/C/3M3 /i3 



and 



where M"^'^ is the charged lepton mass matrix. In a good approximation it 
can be treated as a diagonal matrix M^'^ = Diag{nip} with mf"* being the 
physical lepton masses. Also, one can safely neglect matrix E' compared to 
the other entries of the full mass matrix ( [A.12| ) taking into account smallness 
of the lepton masses. 

Rotation to the mass eigenstate basis 

^(±)^ = (A.16) 
casts the mass matrix in Eq. ( [A.12|) to a diagonal form 

(/\-)* M± (A+)^ = Diag{ml\m^±}, (A.17) 

where mf and m^± are the physical charged lepton and chargino masses. 



2 Appendix B 

Here we give a short account on the results of the approximate diagonalization 
method used in our analysis. 
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2.1 Neutral fermion mixing matrix. 

To leading order in the small expansion parameters ^ defined below, an ap- 
proximate form of the neutral fermion 7x7 mixing matrix introduced in Eqs. 
®, O is 



(B.l) 



Here 



/ii ^2 (Ml + tan^ ^lyMs) sin/5 cosO^Mz . ,^ 
2 UetM^ 

. ^2 (Ml + tan2 ewM2) cos/5 cos^^i/Mz , ,^ 

= YdmT^ ^^-^^ 



with i = 1,2, 3. The determinant of the MSSM neutralino mass matrix ( [A .61) 
is 

r2 ../S/f I -i-„„2 /) n/r t\ /r ..2 



DetM;^ = sin2/?M^/i(Mi + tan^ ^vkM2) - MiMa/i . (B.5) 
>< 4 

diagonal form 



The 4x4 matrix rotates the MSSM neutralino mass matrix to the 



N*M^N^ = Diag{m^^} (B.6) 

where rriy.. are the physical neutralino masses. Thus, to leading order in the 
mixing within the neutralino sector is described as in the MSSM by 

Xk = NknX'n (B.7) 

with Xn = ("^-^'j ""^-^S! -^1 ) -^2) being the weak basis. 

The 3x3 matrix V^'^^ rotates the RPM induced effective neutrino mass 
matrix to the diagonal form 

Vj rueff K = Dtag{0, 0, m,} , (B.8) 

The tree level expression for this mass matrix can be found in Ref. [1^. The 
only non-zero neutrino mass is given by 

, Mi + tan^gH/M2 ,;^,2 
= 4 DetM, • ^^-^^ 
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where 



Let us show an exphcit form of the neutrino mixing matrix 

- sin 6*13 



/ cos 6i3 — sin 613 

sin 623 sin 6^13 cos 623 sin 6^23 cos 613 
sin 6*13 sin 6*23 cos 6*13 cos 6*23 



V 



(B.IO) 



(B.ll) 



where the mixing angles are expressed through the vector A as follows: 

, n Ai , ^ A2 

tan6'i3 = 

VA2 + A3 

described by 



(B.12) 



0\fTAi A3 
The mixing within the neutrino sector to leading order in ^ is 

= V!,tu'^ (B.13) 
with v'^ = (z/e, z/^, Ur) being the weak basis. 



2.2 Charged fermion mixing matrix. 

To leading order in the small expansion parameters and C,^ defined below, 
an approximate form of the charged fermion 5x5 mixing matrix introduced 
in Eqs. ( [A. 161 ) reads 



Vl{1 - 



and 



Here 



92 



V2 DetM,^ 
with 2 = 1,2,3 and 



L* _ fii 92 sin/? cosO^yMz 



/i DetM^± 



A,. 



(B.14) 



(B.15) 



(B.16) 



(B.17) 
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This matrix is much smaller than C,^ by the factor mi/MsusY, where mi and 
MsusY are the lepton masses and the typical SUSY breaking scale Msusy ~100GeV. 
Thus the mixing between (e^^, /ij, r^) and (— 'jA+, H2) described by the off 
diagonal blocks of the A"*" in Eq. (p.l5| ) is small and, therefore, neglected in 
our analysis. 

In Eqs. ( p.l4|) - ([B.15| ) the determinant of the MSSM chargino mass matrix 

is 

DetM^± = Ma/i - sin2/3M^. (B.18) 

The other matrices are defined as follows: 

U*M^±V^ = Dtag{m^±}, (B.19) 
VLM^'^V^ = Dtag{mi^}, 

with M^± and M^') are the MSSM chargino and charged leptons mass matrices 
defined in Appendix A while m^± and m;. are the physical chargino and the 
charged lepton masses. 
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Figure Captions 

Fig.l.: Feynman graphs contributing to neutrinoless double beta (OuPP) de- 
cay for the case of the bihnear R-parity violation, (a) the Majorana neutrino or 
neutralino exchange with two accompanying W-bosons; (b) the gluino-squark- 
squark exchange; (c,d) the neutrino-squark/slepton exchange with one accom- 
panying W-boson. 
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